arXiv:1502.04454vl [physics.optics] 16 Feb 2015 


AXIAL DEPENDENCE OF OPTICAL WEAK 
MEASUREMENTS IN THE CRITICAL REGION 
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Abstract. The interference between optical beams 
of different polarizations plays a fundamental role 
in reproducing the optical analog of the electron 
spin weak measurement. The extraordinary point 
in optical weak measurements is represented by 
the possibility to estimate with great accuracy the 
Goos-Hdnchen (GH) shift by measuring the dis¬ 
tance between the peak of the outgoing beams for 
two opposite rotation angles of the polarizers lo¬ 
cated before and after the dielectric block. Start¬ 
ing from the numerical calculation of the GH shift, 
which clearly shows a frequency crossover for inci¬ 
dence near to the critical angle, we present a de¬ 
tailed study of the interference between s and p po¬ 
larized waves in the critical region. This allows to 
determine in which conditions it is possible to avoid 
axial deformations and reproduce the GH curves. 
In view of a possible experimental implementation, 
we give the expected weak measurement curves for 
gaussian lasers of different beam waist sizes propa¬ 
gating through borosilicate (BK7) and fused silica 
dielectric blocks. 
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I. INTRODUCTION 


The Goos-Hanchen [GH] shift [IH3] surely represents one of the most intriguing research subjects 
appeared in literature in the last decades ima. This shift, which is probably one of the clearest 
manifestations of the evanescent nature of light, represents an additional contribution to the geomet¬ 
rical optical path predicted by the Snell law This quantum effect is still subject of careful and 

broad investigation and continues to stimulate new discussions [T51122j . Of particular interest tothe 
study presented in this paper, it is the GH shift frequency crossover [53]. For incidence angles, 0o, far 
from the critical angle, 9c, it is well know that the GH shift is proportional to the wavelength, A, of the 
optical beam mMsm- For incidence at critical angle the GH shift is amplified by a factor y^Wg/X, 
where Wq is the beam waist. This amplification has been recently obtained analytically by using the 
stationary phase method [251126] and then confirmed by numerical calculations |23j . This frequency 
crossover will play a fundamental role in deriving the expected experimental curves for optical weak 
measurements in the critical (angle) region. 

In a recent interesting experimental paper ED, by using the optical analog [551150] of the elec¬ 
tron spin weak measurement ED, the behavior of the GH shift curve has been reproduced in the 
region in which the incidence angles are far enough from the critical angle to permit some important 
approximations. 

It is important to observe that in the optical analog of the electron spin weak measurement, 
polarized light plays the role of the spin i particles and the laser beam replaces the coherent electron 
beam. Due to the fact that the displacement produced by the optical system is a lateral shift rather 
than an angular deflection, as happens in presence of the Stern-Gerlach magnet ED, we have to 
consider spatial distributions instead of momentum distributions. Notwithstanding the physics is far 
from being the same, by using the needed attention, an optical version of the electron spin weak 
measurement experiment can be constructed |28j . 

A unified linear algebra approach to dielectric reflection, recently appeared in litterature [551130] , 
bases the analogy between weak values and optical beam shifts of polarized waves on the expectation 
value of the Artman operator. Such an operator is shown to be Hermitian for total internal reflection 
and non-Hermitian in the critical region [29]. For the mathematical details, we refer the reader to 
ref. |30j . In our approach, we discuss the optical analog of the electron spin weak measurement, by 
analyzing, as done theoretically in ref. [28] and experimentally in ref. ED, the distance between the 
peaks of the outgoing optical beam. We recall, that in the critical region, due to the breaking of 
symmetry [55], peak and mean value does not necessarily coincide. In view of these comments, our 
discussion can be seen as a complementary work to that one which appears in refs. [291130] . 

To make this introduction and the objective of our analysis clearer to the reader, we recall that 
in the optical analog of the quantum weak measurement [271128] . the parameters which characterize 
the behavior of the peaks distance in the experimental curves are 

€ = Co + Ae = cos(a — /3)/ cos(q! + 13) , (1) 

where a and fd = a + ^ +7o + ^7 are the polarization angles of the first and second polarizer, see 
Fig. 1, and 

^Vgu = Von - Vgu , ( 2 ) 

[s,p] 

where are respectively the GH shifts for s and p polarization. For small rotation angles, i.e. 
Ay <C 1, and for an incoming beam with an equal mixture of polarizations, i.e. a = 7r/4, we have 

Ay 

€o + Ae = tan(yo -|- Ay) Ri tanyo H-^— . 

cos yo 

For incidence angles far from the critical angle, the condition 

Ae > Aycjj/w(z) r; X/w{z) , (3) 


where w( 2 :) = Wq y 1 -f (A x/ttw^) , is easily satisfied and, as we shall see in detail later, the distance 
between the peak of the outgoing beams for two opposite rotations in the second polarizer, i.e. 


1 







= ITT + 7o ± IA 7 I, is given by 

Ab;,,, ~ Ay^^/\Ae\ . (4) 

Consequently, for polarizer rotations which satisfy the constraint ([3]) , the experimental curve of AY^^^ 
reproduces the GH shift curve amplified by the factor l/|Ae|. The GH behavior and its amplihcation 
(far from the critical region) has been recently confirmed in the experimental investigation presented 
in ref. [27] . 

As observed in the begin of this introduction, the frequency crossover in the critical region [23] 
leads to 

A?/gh \Aw(^ > a . (5) 

This critical GH shift behavior stimulates the investigation of what happens for incidence angles 
near to the critical angle, where, due to the amplification y/w{z)/X, the condition ([3]) could be no 
longer valid. This should modify the shape of the experimental curves and a new formula should be 
introduced to estimate the GH shift by the measurement of the peaks distance, AY^^^. In view of a 
possible experimental implementation of optical weak measurements for incidence near to the critical 
angle, we analyze the expected experimental curves for mixed polarized laser gaussian beams, with 
A = 633 nm and Wo = 200, 300, 500 ym, propagating through BK7 and fused silica dielectric blocks. 

The paper is organized as follows. In section H, we give the transmission coefficient for the beam 
propagating through the dielectric structure of Fig. 2 and calculate the axial dependence of the GH 
shift for BK7 (Fig. 3) and fused silica (Fig. 4) blocks. In section HI, we introduce the idea of weak 
measurement in optics and analyze the effect that polarizer and analyzer play on the s and p polarized 
waves of the outgoing beam. For critical incidence, new parameters have to be introduced (Fig. 5). The 
analysis of the distance between the main peak of the outgoing beams for two opposite rotation angles 
of the second polarizer is presented in section IV. There a new analytical relation between the GH shift 
and the peaks distance is also introduced In this section, we also present the expected experimental 
curves for incoming gaussian beams with different beam waists (wo = 200, 300, 500 fim) propagating 
through BK7 (Fig. 6 ) and fused silica (Fig. 7) dielectric blocks. The axial dependence of optical 
weak measurements is clear in the plots and represents one of the important results of our analysis. 
Gonclusions and outlooks are drawn in the final section. 


II. THE GH SHIFT FOR BK7 AND FUSED SILICA BLOCKS 


In order to obtain the mathematical expression for the transmitted beam, , propagating in the 
y-z plane through the dielectric block (see Fig. 1 and Fig. 2), let us first introduce the gaussian wave 
number distribution which determines the shape of the incoming beam, , 


9(0-do) 


kWg 



{kwo f{0 - Qf) j 4 


( 6 ) 


In the electric amplitude expressions the superscript notation distinguishes between s and p polarized 
light. By using the paraxial approximation (fcwo > 10), the incoming electric field, which moves from 
the source laser to the left side of the dielectric block, can be represented by [T31fT4] 




7r/2 

E^E'^^ / de 9 ( 6 - 60 ) exp 
J- 7r/2 


i [6 - 60 ) ky - i k z 


Eo 


exp 


1 + 2G 




w^ + 2z- 


(7) 


For gaussian lasers with a small beam waist with respect to the dimensions of the dielectric block, we 
can use the step technique of quantum mechanics [33H33] and give the Fresnel coefficients in terms of 
the angle 6 , 'ip, and ip (see Fig. 2), 


sin 6 = n sin ip 


and 


‘7’ = 4 ■ 
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The transmission coefficient which characterizes the outgoing beam is obtained by the transmission 
through the left/right sides and the reflection between the up/down sides of the dielectric block. After 
simple algebraic manipulations (for more details see ref. [37)1. we find 


and 


where 


^ncosipcosO j n cos (p — \J\ — n sin^ ip 

^ (^) “ -2" I - 

(cos 6 + n cos tjj) Xn cos ip + — sin^ ip 

4ncos'0cosd ( cos (p — n\/l — sin^ (p ' 


exp[f 


r“ ( 0 ) = 
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exp[i(^3 
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For nsini^ < 1, the outgoing beam, 


E'::\y,z)=E,e 


ik z 


/ dd T (d) g(d - do) exp 

i {9 - 9o) ky - i kz 
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is centered at 

2/snell 


d(t>s 


kd 0 


= cos do 


_ j[£) 

(tan '00 + 1) AB + (tan tpo — tan do) 


It represents the well-known geometrical shift predicted by the Snell law in ray optics. 
For n sin > 1 an additional phase comes from the internal reflection coefficients in 




r r 

[C.C] 


arctan 


^yTAsin^lpA^ 
n cos p 


, arctan 


sin^ p — \ 
cos p 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


and (O, 

( 12 ) 


and a new shift (the well known GH shift) has to be considered. The numerical data for the prop¬ 
agation through BK7 (n = 1.515) and fused silica (n = 1.457) are respectively plotted in Fig. 3 and 
Fig. 4. The data clearly show the amplification for incidence in the critical region and they are in 
excellent agreement with the analytical prediction for incidence far from the critical angle, 
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n sm po — cos^ po 


(13) 


where the Wo dependence disappears [U131 [221 m]. The plots of the GH shift for BK7 and fused 
silica clearly show an axial dependence. This axial dependence, which has been recently investigated 
and recognized as a possible source of angular deviations in the optical path predicted by the Snell 
law [37], has to be seen in the optical weak measurement curves as well. Understanding how this axial 
dependence modifies the optical weak measurement curves in the critical region is one of the main 
objectives of our study. 


III. WEAK MEASUREMENTS IN OPTICAL EXPERIMENTS 


On the basis of the results presented in the previous section, we can approximate the s and p polarized 
outgoing beams as follows 
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(14) 
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where for , which represents the only entry for which we have not a full analytical expression, we 
have to use the numerical data plotted in Fig. 3 and Fig. 4. 

The intensity of the outgoing beam coming out from the dielectric block and moving towards the 
analyzer < z < Zj^ in Fig. 1) is given by 


4 ut (2/> ^out <z<Zj,) = sin a {y,z) + cos a {y, z) 


oc 


r tana exp 


y ^Snell ^GH 

w(z) 


i A(j)(~ 


exp 


[p] ' 

y ^Snell ^GH 

w(z) 


, (15) 


where t = Tg '/ 'Tg'*' and o ~ '(’gh o ■ "I'l*® dependence of r and jj is plotted in 

Fig. 5a (BK7) and Fig. 5b (fused silica). After removing the phase difference between the s and p 
polarized light by the analyzer located at z = and combining s and p polarization by the second 
polarizer located at z = z^, the outgoing beam intensity becomes 
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Observing that 


[O] [p] 

Y = y - 2/s„en - 2 ^ ■ 


tan a tan (3 = 


e - 1 eg - 1 


Ae 


e + 1 Cg + 1 (Cg + 1 ) 

the choice of an appropriate rotation 7 g permits to fix the parameter eg(= tanyg) to 

T - 1 Cg - 1 1 


Eo = 


1 + r 


Cg + 1 


(16) 


(17) 


The angular dependence of 7 g is plotted in Fig. 5b (BK7) and Fig. 5d (fused silica) for incidence angle 
in the critical region. This choice allows to rewritten the outgoing intensity in terms of the parameter 
T and Ae as follows 


^out (Y, Ae) oc 
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exp 
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exp 
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exp 
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(18) 


Finally, noting that, in the critical region, (1 + r) « 4t (see Figs. 5b-d), we can get a further 
simplification of the outgoing beam intensity, 
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(19) 
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IV. THE PEAKS BEHAVIOR IN THE CRITICAL ANGLE REGION 


The starting point in optical weak measurement experiments is to set the angles of the first and second 
polarizers to 

{ ao ) /3o } = { f T + 7o } • 

For this choice (Ae = 0) we find that the outgoing intensity, 


Lut(VO) oc r" exp 


2y' 

w^(z) 


( 20 ) 


is a symmetric function with two peaks centered at = ± w(z)/-\/2 and a minimum centered at 
Lmin = 0. By changing the angle of the second polarizer from /3o to /3o + Ay, we break the symmetry. 
In terms of Ae = Ay / cos^ yo, we find 


V>in(Ae) = 


Ae 

^Vgh 


w=(z) 


( 21 ) 


and 


y (Ae) = 

max ^ ' 


- Ae ± +2[Ayl^/ w=(z)] 

2AyGH 


w=(z) . 


( 22 ) 


It is clear that for positive Ae (anti-clockwise rotation Ay around yo), the y^^^(|Ae|) represents the 
position of the main peak of the outgoing beam. For negative Ae, the main peak is instead centered 
at |Ae|). By using Eq. (1^ . the distance between these peaks is given by 


Ay =y (|Ae|)-y (-|Ae|) = 

max max '-I I / max k I I / 


- |Ae| + ^\Ae\^ +2[Ayl^/ w=( 2 ;)] 


Aj/c 


\z) . (23) 


In the region 0 < |Ae| < Ay^^/w{z), we find 


V2w{z) < AY^^^ < (\/3-1)w(2;) . 


(24) 


This clearly shows that by increasing the value of |Ae|, we reduce the distance between the peaks. 
For |Ae| > Ay^^/w{z), 

Ay_ « Ay^,/\Ae\ . (25) 

For incidence angle far from the critical region, due to the fact that the GH shift is proportional to 
the wavelength of the laser beam, the condition 


A^gh 

w(z) 


A 

w{z) 


< |Ae| 


can be easily satisfied. Thus, far from the critical region, the experimental curves of AY^^^ reproduce 
the GH curves amplified by the factor l/|Ae|. 

In the critical region, the frequency crossover and the axial dependence, showed in Fig. 3 (BK7) 
and Fig. 4 (fused silica), play against the validity of the constraint |Ae| :#> Ay^^/w{z). This means 
that in this region, the experimental curves of Ay^^^ do not necessarily reproduce the GH curves. 
The expected experimental curves of the peaks distance are plotted for different value of | Ae| in Fig. 6 
(BK7) and Fig. 7 (fused silica). The plots confirm that, at critical incidence, the amplification does 
not reproduce the l/|Ae| proportionality. The axial dependence is removed by increasing the beam 
waist Wq. 

For experimental use, it is convenient to express the GH shift, Ay^^, in terms of the experimental 
quantity AY^^^. From Eq. (l23l) . we obtain 


Aj/gh 


2|A£|w-(z) 

2w^(z)-Ay" 

' ' max 


(26) 
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The error on the GH shift is given by 



r^(|Ae|)] 
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Recalling that for Ae = 0, the distance between the peaks gives a direct information on the beam 
waist, = •\/2w(z). we can use 




in the previous error formula and obtain 
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In the region w( 2 ;) < < •\/2w(z), we find 
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2 Ay /w(z) 

max' ' ' 


2wRz) + Ay" 

^ ' m 


< oo . 


To avoid great standard deviations, we have to work in the region Ay^^^ < w(z), where 


K|Ae|) 

|Ae| 
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Ay„. 


^ g(AyGH) ^ 


Avc 


K|Ae|) 

|Ae| 


13 


crfAF ) 

Ay 


From the previous condition on Ay^^^, by using Eg. (1261) we obtain 

AycH 


|Ae| > 


2 w(2;) 


(27) 


(28) 


The choice of |Ae^j„| = Ay^^/2w{z) in the second polarizer thus represents an additional experimen¬ 
tal constraint to avoid great standard deviations. 


V. CONCLUSIONS AND OUTLOOKS 


The possibility to use the weak measurement of the electron spin component [31] in optics [28112^ 
have recently stimulated the realization of an experiment m, based on the interference between 
different polarizations, in which the GH shift curves are reproduced in the region of the validity of the 
standard analytic formula (IT^ . Nevertheless, the analytical shift (fT3l) diverges when the incidence 
angle approaches the critical angle. In a recent paper |23j , this divergence was removed and an analytic 
formula, valid for 2 z <C kwl, proposed for the GH shift at critical angle, 


r [«i [p] I 

\/A:w(z) / 

1 Von > ^GH f - 

^ •' cri 

k V 


2^27 


^2 - n" + 2 Vn" - 1 
n" — 1 -I- Vn" — 1 




(29) 


This closed formula, which is in excellent agreement with the numerical data plotted in Fig. 3-c and 
Fig. 4-c, clearly shows the crossover frequency at critical angle. The amplification y^kw{z) at critical 
angle suggested studying with more care the peaks behavior of the beam in optical weak measurements 
for incidence within the critical region. Indeed, in such a region, due to this amplification, the condition 
I Ae| ^ Ay^^/w(z) and the consequent proportionality amongst the experimental curves of the peaks 
distance and the GH curves are no longer valid, see for example Fig. 6-c and Fig. 7-c. 

In our study, we have also found an axial dependence in optical weal measurements. This axial 
dependence can affect the experimental curves and, for small second polarizer rotations and small 
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values of the beam waist, produces a practically flat region, see the plots in Fig. 6-a and Fig. 7-a for 
|Ae| = 0.01. To minimize the axial dependence, we have to work with laser beam with Wq > 500/rm. 
It is important to be observed here that, also for Wq = 500 ^m, the curve amplification l/|Ae| can be 
only reproduced far from the critical region. 

In view of a possible experimental analysis of the study presented in this article, we have also 
estimated in which region we reach the better standard deviation for Aj/^g in the critical region. By 
using the second polarizer angle constraint, Ea. (l28l) . and the analytical formula for the GH shift at 
critical angle, Eo. d^ . we find 


|Ae| > 



a/ 2 - n" + 2 Vn" - 1 

— 1 + 'Jn 


— 1 

\/kw{7j 


(30) 


This implies, for laser beams with Wg = 500/im and camera at z < 50 cm, |Ae^i,j| ~ 0.015 both for 
BK7 and fused silica dielectric blocks. 

We conclude this work, by observing that our analysis does not take into account cumulative 
dissipations and imperfections in the dielectric prism (such as the misalignment of its surfaces) and 
the beam reshaping caused by interference . A phenomenological way to include misalignment effects 
is given in ref. m- An interesting discussion on the origin of negative and positive lateral shifts in a 
dielectric slab is investigated in ref. [38] form the viewpoint of the interference between multiple light 
beams. 
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EXPERIMENTAL SETUP 


CAMERA 



Figure 1: Experimental layout. A schematic representation of the optical weak measurement 
experiment for the observation of the transversal distance between the main peak of the outgoing 
beams for two opposite rotations, iAy, of the second polarizer. The incoming beam, which, after 
the first polarizer [a = 7 r/ 4 ), has an equal mixture of s and p polarized waves, passes through the 
dielectric block < z < and then passing through the analyzer in z = loses the global 
A^!>qh phase). Optical weak measurements are done by changing the rotation angle in the second 
polarizer (/3 = 37 r /4 + 7 q ± | Ayj). The angle 70 is fixed to obtain, for A 7 = 0, an outgoing beam with 
two identical maxima centered at ± w(z)/y^. 


10 



































Figure 2: Geometrical path and Goos-Hanchen shift. Schematic diagram of the dielectric block 
analyzed in this paper. In (a), it is shown the geometrical path predicted by the Snell law, Eq. CD- 
For (/? > (/?c, an additional phase, coming from the Fresnel reflection coefficients at the down and up 
interfaces, has to be considered. This phase is responsible for the addition shift, Eq. (I13L known as 
Goos-Hanchen shift and shown in (b). 
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Figure 3: The GH shift curves for BK7 blocks. The numerical data for the transversal 
shift of laser gaussian beams, passing through a BK7 dielectric block, are plotted, in the axial range 
10cm < z < 15cm for different beam waists Wq = 200/im (a), 300^m (b), and 500/im (c). The 
crossover frequency at critical angle is clear from the plots and the axial dependence represents an 
additional phenomenon to be considered in optical weak measurements. 
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Figure 4: The GH shift curves for fused silica blocks. The numerical data for the transversal 
shift of laser gaussian beams, passing through a fused dielectric block, are plotted, in the axial 
range 10cm < z < 15 cm, for different beam waists Wj, = 200^m (a), 300/rm (b), and 500/rm (c). 
The crossover frequency at critical angle is clear from the plots and the axial dependence represents 
an additional phenomenon to be considered in optical weak measurements. 
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Figure 5: Angular dependence of x, and 70 . The angular dependence of r (ratio between 

the modulus of the amplitudes for s and p polarized light) and (global phase difference between 

s and p polarized waves) are plotted in (a) for BK7 and (c) for fused silica blocks. The fact that, in the 
critical region, 4r/(l + r) is practically equal to one, it is very useful to simplify the expression for 
the outgoing beam, see Eq. (Iig. The numerical data for 70 permits to calculate the second polarizer 
angle (/3o = 37 r /4 + 70 ) for which we find an outgoing beam with two identical maxima centered at 
±w(z). 
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Figure 6: Optical weak measurements curves for BK7 blocks. The expected curves for the 
distance between the main peak of the beams coming out from a BK7 dielectric block and passing 
through the second polarizer for two opposite rotations, |Ae| = IA 7 I/cos^ 7 o, are plotted in the axial 
range 10cm < z < 15 cm, for different beam waists Wj, = 200/rm (a), 300^m (b), and 500/rm (c). 
From the plots, it is clear that to improve the crossover frequency and to reduce the axial dependence, 
we have to work with Wq > 500 /rm. Note that, also working with Wq = 500 /rm, the curve amplification 
l/|Ae|, valid for incidence far from the critical region, is lost when the incidence angle approaches the 
critical one. 
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Figure 7; Optical weak measurements curves for fused silica blocks. The expected curves 
for the distance between the main peak of the beams coming out from a fused silica dielectric block 
and passing through the second polarizer for two opposite rotations, |Ae| = IA 7 I/cos^ 7 o, are plotted 
in the axial range 10cm < 2 < 15 cm, for different beam waists Wq = 200/rm (a), 300p,m (b), and 
500 /im (c). From the plots, it is clear that to improve the crossover frequency and to reduce the axial 
dependence, we have to work with Wq > 500 /rm. Note that, also working with Wq = 500 /rm, the 
curve amplification l/|Ae|, valid for incidence far from the critical region, is lost when the incidence 
angle approaches the critical one. 
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